The Lyapunov spectrum of the many-dimensional dilute random Lorentz gas 
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For a better understanding of the chaotic behavior of systems of many moving particles it is useful to 
look at other systems with many degrees of freedom. An interesting example is the high-dimensional 
Lorentz gas, which, just like a system of moving hard spheres, may be interpreted as a dynamical 
system consisting of a point particle in a high-dimensional phase space, moving among fixed scat- 
terers. In this paper, we calculate the full spectrum of Lyapunov exponents for the dilute random 
Lorentz gas in an arbitrary number of dimensions. We find that the spectrum becomes flatter with 
increasing dimensionality. Furthermore, for fixed collision frequency the separation between the 
largest Lyapunov exponent and the second largest one increases logarithmically with dimensional- 
ity, whereas the separations between Lyapunov exponents of given indices not involving the largest 
one, go to fixed limits. 

PACS numbers: 05.45.Jn 



I. INTRODUCTION 

Many studies have been done on the chaotic properties 
of the Lorentz gas Q, 0, El • It is a simple model which 
consists of a single particle moving freely between elastic 
spherical scatterers. The scatterers can either be placed 
at random or in a lattice configuration. The Lorentz gas 
provides a physical system yet is still simple enough to 
allow for exact calculations of many properties. This sim- 
plicity is partially due to the fact that the Lorentz gas 
contains only one moving particle, and therefore does not 
have many dynamical degrees of freedom. Systems with 
more degrees of freedom, such as moving hard spheres or 
disks, have also been studied frequently. Extensive sim- 
ulation work has been done on their Lyapunov spectrum 
0- HI 01 j and for low densities analytic calculations have 
been done for the largest Lyapunov exponent 
the Kolmogorov-Sinai entrop y H and for the smallest 
positive Lyapunov exponents [ll], E3 ■ 

From the viewpoint of dynamical systems theory the 
Lorentz gas and hard sphere systems are very similar, as 
noted already many years ago by Sinai [l^. Both sys- 
tems may be viewed as "billiards" , i.e. systems consisting 
of fixed obstacles in a, mostly high-dimensional, configu- 
ration space, among which a point particle moves elasti- 
cally. In the case of the Lorentz gas these scatterers are 
(hyper) spheres, in that of the hard-sphere system they 
are hypercylinders. In either case the convexity of the 
scatterers makes the system strongly chaotic. In several 
respects the Lorentz gas is much simpler than the hard 
sphere system. First of all the scatterers for the Lorentz 
gas are invariant under rotations in configuration space, 
which, as we will see, simplifies calculations enormously. 
Further, the uniform convexity of the Lorentz gas scatter- 



ers, in contrast to the hypercylinders of the hard-sphere 
systems, strongly simplifies proofs of ergodic and chaotic 
properties [T^J. Yet we think it is of interest to perform 
an explicit calculation of the full Lyapunov spectrum of 
a high-dimensional dilute Lorentz gas. It is interesting 
to see the similarities as well as the differences between 
the Lorentz gas and the hard-sphere spectra. In addition 
the methods used here may well be amenable to refine- 
ments, so as to make them applicable to systems of many 
moving particles. 

In this paper, we study the behavior of a dilute, ran- 
dom, non-overlapping, Lorentz gas in an arbitrary num- 
ber of dimensions, d. For large d this system has many 
degrees of freedom, while, largely due to the spherical 
symmetry of the scatterers, it is still possible to do exact 
calculations. Here we calculate the full Lyapunov spec- 
trum in the absence of any external fields. 

As a preparation we introduce Lyapunov exponents 
and the Lorentz gas in section ITT1 and we discuss the low 
density approximation. In section ITTT1 we derive an inte- 
gral expression for the spectrum of Lyapunov exponents. 
Then, in section ITVl we approximate this expression for 
high-dimensional systems and investigate its properties. 
In section[V]we discuss the results and make comparisons 
to the hard-disk Lyapunov spectrum. 



II. LYAPUNOV EXPONENTS AND THE 
LORENTZ GAS 

Consider a system with an A/"-dimensional phase space 
r. At time t — the system is at an initial point j Q 
in this space. It evolves with time, according to 7(70, t). 
If the initial conditions are perturbed infmitesimally, by 
($70, the system evolves along an infmitesimally different 
path 7 + <57, specified by 
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-57(70, t) = M 70 (t) -Sj 



(1) 



in which the matix M 7o (t) is defined by 



,(*) 



dj 



(2) 



The Lyapunov exponents are the possible average rates 
of growth of such perturbations, i.e., 



K = lim 7 log \m(t)\ , 



(3) 



where is the i-th eigenvalue of M T0 (i). If the system 
is ergodic, it comes arbitrarily close to any point in phase 
space for all initial conditions except for a set of measure 
zero. Therefore, the Lyapunov exponents are the same 
for almost all initial conditions. We will order the ex- 
ponents according to size, with Ai being the largest and 
Xj^f the smallest, as is the convention. For each exponent 
there is a corresponding eigenvector of M 7o (t). 

The dynamics of a purely Hamiltonian system are com- 
pletely invariant under time reversal. Also, for ergodic 
Hamiltonian systems, due to the incompressibility of flow 
in phase space, the phase-space attractor is invariant un- 
der time reversal. Therefore, every tangent-space eigen- 
vector that grows exponentially under time evolution, 
shrinks exponentially under backward time evolution. As 
a consequence, since the Lyapunov spectrum does not 
change under time reversal, there is a negative exponent 
of equal absolute value for every positive Lyapunov ex- 
ponent. This is called the conjugate pairing rule. 



A. The Lorentz gas 

We will consider the dilute random Lorentz gas. It is a 
system with a fixed number N of randomly placed spher- 
ical scatterers with diameter a at a (small) density n in 
d dimensions. To be specific, we will restrict ourselves 
to the case of non-overlapping scatterers, where, a priori, 
each configuration without overlap among any two scat- 
terers is equally likely. We will assume that the system 
is large. As long as it is finite our preferred boundary 
conditions are periodic ones, but our considerations al- 
low taking the infinite system limit, at fixed n, without 
any problem. A single point particle moves between the 
scatterers, with speed v, undergoing a specular reflection 
at each collision. For given scatterer positions, the phase 
space is represented by the position and velocity of the 
point particle, 7 = (r, v). The tangent phase space at 
any point in phase space can be represented by the per- 
turbations in these quantities, £7 = (<5r,<5v). 

Between collisions with the scatterers, the particle 
moves freely, so r = v and v = 0. At a collision with col- 
lision normal <r, the particle is reflected by a scatterer, as 
shown in Fig.^ with resulting velocity v' = (1 — 2aa) - v. 
Here 1 is the d x d identity matrix. From Eq. the tem- 
poral behavior of $7 can be derived. During free flight, 



St 
Sir 



8v 
. 



(4) 
(5) 




FIG. 1: Geometry of a collision. The collision normal a is the 
unit vector pointing from the center of the scatterer to the 
moving point particle. 



At a collision, Sv is also reflected. On the tangent trajec- 
tory, the perturbation of the precollisional position with 
respect to the reference trajectory, leads to a perturba- 
tion in the collision normal and the collision time, which 
shows up in the perturbation of the postcollisional ve- 
locity. The perturbation of the collision time also leads 
to a slight deviation of the postcollisional position differ- 
ence (measured at the instant of collision on the tangent 
trajectory) from the precollisional one (measured at the 
instant of collision on the reference trajectory). As a 
result, the tangent phase space vectors are transformed 
according to 



Sr' = (1 - 2aa) ■ St . 
Sv' = (1 - 2aa) • Sv 



2Q-Sr 



(6) 
(7) 



in which the collision matrix Q is defined by 

If v 2 \ 

Q = - V(j - (tv H r <7<7 - (v • (j)l ■ (8) 

a \ v • a J 

From these equations it follows that, if Sr and Sv are both 
parallel to v before the collision, Sr' and Sv' are parallel 
to v' after the collision and their absolute values are the 
same as before. These are two linearly independent per- 
turbations, giving rise to two zero Lyapunov exponents. 
They result from time translation invariance and from in- 
variance of the trajectories in configuration space under 
a scaling of the velocity. All other Lyapunov exponents 
for the Lorentz gas are non-zero. As a consequence of 
the conjugate pairing rule d — 1 of them are positive and 
the remaining d — 1 are negative with the same absolute 
values. 

B. The largest Lyapunov exponent at low densities 

To calculate the non-zero Lyapunov exponents, one 
needs to consider what happens to an initial perturba- 
tion (57(0) in tangent space in the limit of infinite time. 
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As an introduction we first review the calculation of the 
largest Lyapunov exponent at low scatterer densities. 

An initial perturbation which is not parallel to v gener- 
ically has a non-vanishing component along the most 
rapidly growing eigenvector of the time evolution oper- 
ator in tangent space. Therefore its evolution for long 
times will be dominated by the largest Lyapunov expo- 
nent. To calculate this time evolution it suffices to con- 
sider the growth of the projection of the growing vector 
onto a subspace of tangent space. It turns out convenient 
using the projection onto <5v for this. 

Define Sri and Svi as the tangent space vectors just 
after collision i, with collision normal a"j, occurring a time 
n after collision i — 1. Let 9i be defined by 

(9) 



-wcosf 



Before collision i one has 



5r 7 = Sr l _ 1 +t z Svi-i, 



(10) 



where Sr~ is used to indicate the perturbation in the 
position just before collision i. From Eqs. J7J) and (JSJ it 
follows that, typically, after the (i — l)-th collision <5vj_i 
is of the order of Sri_iv/a. At low densities the mean 
free time is of the order of 1/ (nav). Therefore, to leading 
order in the density the first term on the right hand side 
of Eq. Ijl 0(1 may be neglected. Similarly, in Eq. Q the 
first term on the right hand side becomes negligible at 
low density, and Eqs. ((10(1 and (JJJ may be combined into 



5vi = 2nQi ■ 5vi 



(11) 



The contributions to the Lyapunov exponents of the 
terms neglected in this approximation are at least one 
order of n higher than the terms of leading order pj . In 
addition, in this approximation, the time reversal symme- 
try has been destroyed, hence only the positive Lyapunov 
exponents can be calculated. However, for the limit of 
density going to zero, the results for these will be exact. 

The action of Qi on #Vj may be described in the fol- 
lowing way. Working on the component along Vj it yields 
zero. It multiplies the component normal to v in the 
plane through Vj and a with a factor v /(a cos 9i) and 
rotates it to the direction in this plane normal to v',. 
Finally it multiplies all other components of <5vj with 
v cos 9i/ a. Define the unit vector orthogonal to v in the 
plane spanned by v and a as 

(1 - VjVi) ■ frj 

I sin 9i\ 

One may rewrite Qi as 

Qi = - ( cos#j(l - isiisi - piPi) H ^—p-p'iPi J -(13) 

a \ cost', / 

Combining this with Eqs. ((6I7() one finds that the velocity 
deviations to leading order in n evolve according to 



6vi = 



2vn 



cos0j(l - v,Vj - pi^ 
1 A 



cos 9 



■Pi Pi 



■ Sv 



i-l 



(14) 



The largest Lyapunov exponent may now be calculated 
to leading order in the density, as 



\i = Vd( log 



|<5v, 



(15) 



where Vd is the average collision frequency for the system 
and the brackets indicate an average over the collision se- 
quence, which will be discussed in more detail in Sec. lIIII 
In 2 and 3 dimensions these calculations have been done 
before [lj. The result in d dimensions will appear as a 
special case of the calculations presented in the next sec- 
tion. 



III. PARTIAL STRETCHING FACTORS 

In standard terminology, the stretching factor is de- 
fined as the factor by which the expanding part of tangent 
space expands over a time t. This quantity can be used 
to calculate the Ruelle pressure as well as the sum of the 
positive Lyapunov exponents, equaling the Kolmogorov- 
Sinai entropy in systems without escape 0, • 

We may define the partial stretching factor As(r, v, t) 
of a p-dimensional subspace S of the 2d-dimensional tan- 
gent phase space as the factor by which the volume of an 
infinitesimal p-dimensional hypercube in this subspace 
has increased after a time t. Unless S is orthogonal to 
some eigenvector associated with one of the p largest Lya- 
punov exponents, the partial stretching factor for very 
long times will be dominated by the p most unstable di- 
rections in tangent phase space, in other words, by the 
p largest Lyapunov exponents. Explicitly, one has the 
identity 



1 



p 

V Xi = lim - log A s (r, v, t) . 

' — ' t— >oo t 

1=1 



(16) 



As in the case of the largest Lyapunov exponent, where 
we could consider the long time growth of a basically ar- 
bitrary vector in tangent space, we may choose the sub- 
space S in the way that is most convenient to us. And, as 
before, we choose S as a subspace of the space spanned 
by velocity deviations perpendicular to v. 

The partial stretching factor just after collision N is 
the product of the partial stretching factors of the colli- 
sions 1 through N. These depend on the relative orien- 
tations of v, <7, and the image of S. One can write 



JV 



A s (r,v,t N ) = Y[A$H v ,n,0i,(Xi) 



(17) 



in which a, is the projection angle of pi onto the im- 
age of S after the (i — l)-th collision. The subspace S 
can be split into a (p — l)-dimensional subspace normal 
to pi and a 1-dimensional subspace spanned by the pro- 
jection of pi onto S. From Eq. 1(14(1 one finds that the 
former contributes a factor of (r, cos 9i) p ~ 1 to the partial 
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stretching factor. The projection of pi onto the image of 
S can be split into components perpendicular and par- 
allel to pi. The former grows with 2vn cos Qija and the 
latter with 2vTi/(acos8i). The partial stretching factor 
thus becomes 



AV\v,Ti,6i,ati) 



cos p_1 9i 



x\ (sin cti cos Oi)" 1 



cos Cti 
cos 9i 



(18) 



From this expression one can calculate the Lyapunov ex- 
ponents and obtain asymptotic approximations for high 
dimensionality. 



since the probability of a return to a scatterer decreases 
rapidly. In this approximation the time of free flight be- 
tween consecutive collisions is distributed exponentially. 
In addition the distribution of the angles 9 and a is inde- 
pendent of that of the free flight time, the direction of the 
incident velocity and the orientation of the precollisional 
image of S. The probability density for colliding at angle 
9 is proportional to the differential cross section. In the 
free flight time distribution it contributes a factor of 



x v O d -ism d - 2 9 cos 9 d9dt 



(21) 



Here O m is the (m — l)-dimensional surface area of the 
77i-dimensional unit sphere, i. e., 



Or, 



2tt- 



(22) 



A. Lyapunov exponents 

The sum of the p largest Lyapunov exponents can be 
calculated by substituting Eq. I|17(l into Eq. (11611 . 



p 



i=l 



N 



lim J- VlogAW^,^,^,^) , (19) 

jv— >oo Ijy * — ' 
i— 1 



where tpj is the time at which the iV-th collision occurs. 

In the low density limit, the collisions in the Lorentz 
gas are uncorrelated and therefore the time average in 
Eq. i|19|) can be replaced with an ensemble average, 



Finally, the probability distribution of the projection an- 
gle a of p onto the image of S after i — 1 collisions may 
be identified with the fraction of the (d — l)-dimensional 
unit sphere that has a projection angle between a and 
a + da, viz. 



O d . 



1—pOp • d — 2 — Ti 

— — - sin p 



p(a) = ^ O d -i 
5(a) 



a cos p 1 a if p < d — 1 
]£p = d-l 



(23) 



where 5(a) is the Dirac delta function. Combining 
Eqs. (|21|l and 11': ill with the exponential distribution of 
the free flight times, one finds, for p < d — 1, 



5> = 



dr d9 da 
(t Jo Jo 



x v p (r,e,a;v)logA^(v,T,6,a) . (20) 



Here v p (t, 9,a;v) is a probability distribution, describ- 
ing the probability density per units of time and angle 
for collisions with the parameters r, 9, and a at a given 
velocity v. The index p is attached to remind of the de- 
pendence of the distribution of a on the dimensionality 
of the subset 5*. A more rigorous proof of the theorem 
used to derive Eq. I|2U|) can be found in Crisanti, Paladin, 
and Vulpiani 



m 



f p (r,6,a;v) = 

2na d ~ 1 v Od-i-p O p v d exp -v d r 
x sin^ 2 9 cos (9 sin d - 2 " p a cos 25 " 1 a , (24) 

where the average collision frequency has the explicit 
form 



na vO d - 
d- 1 



(25) 



Note that O d -\j(d — 1) is equal to the volume of the 
(d — l)-dimensional unit sphere, so a d ~ x O d -\l (d — 1) is 
the total cross-section for a collision with a scatterer. 



B. The distribution of free flight times 

IV. THE SPECTRUM 

To lowest order in the density of scatterers collisions 
are uncorrelated; effects of recollisions appear in the Lya- Substituting Eqs. (fTH|l and ij^ into Eq. (|2Tj|) yields an 

punov exponents only at higher orders. With increasing expression for the sum of the p largest Lyapunov expo- 
d these higher density corrections even become smaller, nents, 



5 



5> 



c?t / dfina v Od-\Vd exp(— v&t) sin "0cos( 

I! ./O 



/ 2vr 

p log ( ) + (p — 1) log cos 6> 



dft / dana d " 1 ?;O d _i_pOpsin d " 2 ecos^siii d_2 " p acos ? '" 1 alog 

I 



. . m2 /cos a 

(sin a cos ft) + - 

V cos 9 



(26) 



The first integral can easily be carried out analytically. 
The second integral is more difficult. It can be simplified 
in special cases, as described in subsection IIV Al or ap- 
proximated, as described in the remainder of this section. 



A. The KS entropy 



tangent space vector is realized for every collision. The 
minimum growth is the growth in the directions perpen- 
dicular to a and v. This yields, 



Ad-i > A_ , 



(29) 



By taking p = d— 1, one can calculate the sum of all the 
positive exponents, the Kolmogorov-Sinai (KS) entropy. 
Since the distribution of a now is a delta function, the 
second term in Eg. 1261 becomes 



dOna vOd-i-pO p sin 9 cos ( 
sin d_2_p a cos p_1 a log cost 



(27) 



Performing the integrals yields, 
d-i 



»=i 



avd 
2v 



x{ 2(d- 1) 



(d-3) 



-log 



(0) 



f na d O d ^ \ 
\2(d-l)J 



+ 7 



(28) 



Here t/j^(x) is the (i + l)-th derivative of the logT(a;); 
ij)^{x) = [logr(a;)]( i+1 ). This reproduces the results of 
Van Beijeren, Latz, and Dorfman 



B. Lower bound 

A lower bound for the Lyapunov exponents can be 
derived by assuming that the minimum growth of the 



with A_ defined by 



2 log 



V 2v 



37 — ip 



(0) 



(30) 



The function ip^(d + 1/2) for large d behaves as logd. 
Therefore the log v d term dominates the behavior of A_ , 
and A_ for large d behaves asymptotically as i>d\ogVd- 

In a system of high dimensionality, the exponents are 
dominated by the directions perpendicular to v and a. 
They should behave as A_ with a small correction. 



C. High dimensionality 



The integrals in Eq. (|26|l can be estimated for large 
d and arbitrary p. For large numbers of dimensions the 
distribution of 9 is sharply peaked near 9 = tt/2. The 
argument of the logarithm in the second term of Eq. I|26|) 
is therefore dominated by the 1 / cos 2 9 term. The range 
of a where the other term dominates is of the order of 
1/d, and therefore that term can be neglected. This leads 
to 



2jAi w p\- 



i=l 



d0 / dana d - 1 vO d -i-pO p sm d - 2 9 cosft sin d ~ 2 ~ p a cos^ 1 a (-2 log cos 9 
Jo 



which yields, after performing the integrals, 
p 



5> 



PA_ + T 



27-1/" 



(0) 



d-l 



2i> 



(o) 



d + 1 



* m (I) 



log cos a) j(31) 



(32) 
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For the largest exponent this means that the behavior 
for large d is approximately 



Ai « A_ + 



log d + 7 - | log 2 



(33) 



This behavior is shown in the uppermost curve in Fig. 01 
The smallest exponent is equal to the KS-entropy mi- 
nus the sum of the first d — 2 exponents. In the high 
dimensionality limit the smallest positive exponent be- 
haves as, 



Ad-i ~ A_ 



2d 



(34) 



This is illustrated in Fig. The p-th exponent can be 
calculated by subtracting the expression for the sum of 
the first p — 1 exponents from that for the first p ones. 
This results in 



A P «A_ + ^ 



1> 



(o) 



(f) 



if, 



(0) 



p-1 



For the first few exponents this yields 



A 2 
A 3 
A 4 



A_ 
A_ 
A_ 



^dfog2 , 
P d (l-log2) , 
v d (\og2-\) 



(35) 



(36) 
(37) 
(38) 



As a consequence of this, for fixed collision frequency, 
(Ai — X—)/9d grows logarithmically with dimensionality, 
whereas (A p — A_)/P<2 withp > 1, approaches a limit that 
is independent of d. This too is illustrated in Fig. [3] 



V. DISCUSSION 

The integrals in Eq. H26I) can be performed numeri- 
cally. The results are displayed in Figs. 121 ol along with 
the limiting behaviors for large dimensionality, which we 
discussed already in the previous section. 

The offset of all Lyapunov exponents, A_, depends on 
both density and dimensionality. For large d its mag- 
nitude is determined primarily by dimensionality, unless 
density is so low as to satisfy 



na < 



0,j 



1 



(39) 



It is interesting to compare the Lyapunov spectrum of 
a high dimensional Lorentz gas to that of a system of 
moving hard disks, as computed by Dellago et al. [l7j . 
Fig. shows such a spectrum for 750 particles in 2 di- 
mensions at a density of n = 0.1 /a 2 . One immediately 
notices that the Lorentz gas spectrum is much flatter 
than that of the hard disk system. The explanation for 
this is that for the hard disk system only a few (4 to be 
specific) components of the velocity in phase space are 
changed at each collision, whereas in the Lorentz gas all 
components are involved, though some will be increased 
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FIG. 2: The scaled spectrum of positive Lyapunov expo- 
nents for several values of the dimensionality d at density 
n = 0.01/a 2 . 




FIG. 3: The four largest Lyapunov exponents as functions of 
the dimensionality. Their limiting behavior is indicated with 
lines. The second, third, and fourth exponent converge to 
horizontal lines, as calculated in Eqs. I3bl381 . 



more than others. A greater similarity may be obtained 
by replacing the spherical scatterers in the d-dimensional 
Lorentz gas by randomly oriented hypercylinders. In- 
deed our calculations indicate that the spectrum ob- 
tained for this system resembles much more that of the 
hard disk system, but there remain significant differences 
due to the randomness of the cylinder orientations (for 
hard disk systems the scatterers in phase space are hy- 
percylinders of very specific orientations). In addition, 
due to the absence of velocities for the scatterers, the 
rf-dimensional Lorentz gas will never exhibit a branch of 
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FIG. 4: The discrete derivative of the largest Lyapunov ex- 
ponent with respect to log d. The derivative converges to | , 
as calculated in Eq. (1331 . 




FIG. 5: The smallest positive Lyapunov exponent as a func- 
tion of the dimensionality. 



Goldstone modes with Lyapunov exponents approaching 
zero as the inverse of system size, in contrast to hard 
disk systems 0, 0, However, there are also some 
similarities between the two (positive) spectra: first of 
all both become increasingly flat with increasing index 
(apart from the Goldstone branch for the hard disks). 
Secondly both become very steep near the largest expo- 
nent. It is especially remarkable that, for fixed collision 
frequency, the difference between the largest exponent 
and the next largest one in the Lorentz gas increases 
logarithmically with the number of degrees of freedom, 
whereas all the subsequent differences, in the positive half 
of the spectrum, approach fixed limiting values. Whether 




500 1000 1500 

i 

FIG. 6: The spectrum of Lyapunov exponents from sim- 
ulations 

S S Gj of 750 hard spheres in 2 dimensions at 
density n = 0.1/a 2 in a rectangular box of dimensions 
IQ/y/n x 75/y/n, with periodic boundary conditions. Only 
the positive exponents are plotted, since, by the conjugate 
pairing rule, the negative spectrum is exactly the opposite. 
The kinetic energy per particle has been scaled to 1. 



something like this also happens for hard disks is not 
known at present, though there have been conjectures of 
similar behavior by Searles et al [lflj . 



VI. CONCLUSION 

In this paper we have calculated the full Lyapunov 
spectrum for the dilute random Lorentz gas in an ar- 
bitrary number of dimensions. We have found analyti- 
cal expressions for the behavior of the spectrum for sys- 
tems with many degrees of freedom. The spectrum be- 
comes flatter with increasing dimensionality or decreas- 
ing density. The separation between the largest and sec- 
ond largest exponent, expressed in units of the collision 
frequency, increases logarithmically with dimensionality. 

Perhaps a similar approach may be applied to the Lya- 
punov spectra of systems of many particles, such as freely 
moving hard spheres. However, in the present case we 
could take advantage of the property that the partial 
stretching factor is distributed independently of the sub- 
space being stretched 0. This will not be true any 
more for systems of many moving particles, for which 
the dynamics are not invariant under arbitrary rotations 
in configuration space. Some new ideas will be required 
here. 
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